1. Introduction 
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Phase transitions are often described in terms of an order parameter, which 
in the electroweak theory is the background Higgs field. The states with lowest 
energy are conveniently analyzed with the help of the effective potential, which 
^T) ". has lately been the subject of intensive studies in connection with baryogenesis in 
the electroweak theory [1]. It is however rather obvious that to investigate the 
dynamics of the order parameter in detail one should use the complete effective 
action, which includes also terms depending on the derivatives of the background 
field. Although the ground state has a constant value of the background field, 
locally there are fluctuations about this value, and they might turn out to be 
important for local phenomena such as bubble nucleation in first order phase 
^ ! transitions. 

To study the evolution of background fields one should solve the equation 

O 

of motion for the fluctuating background as derived from the effective action, 

Q\ \ 

including possible dissipation. As the full effective action is an object that is very 
difficult to estimate, one might resort to a derivative expansion of the effective 

P+ 

action, cut it off after the first few terms and solve the resulting equation of 



a: 
p. 



motion. This approach does not always work at finite temperature because the 



^ derivative expansion does not necessarily exist [2] in the neighborhood of the origin 

of the (po, p)-plane. 

If the background field is assumed to be close to a constant value so that we 
may consider just small space-time dependent perturbations about that constant 
field, it is possible to expand the effective action in terms of the perturbation 
while retaining all the derivatives. In this approach one is able to see not only 
the dispersion of the fluctuation but also which fields, or modes of the field, 
are unstable and what are their actual decay rates. In what follows we shall 
adopt this approach for A</> 4 -theories. We shall compute the effective action 
in this approximation up to two loops and study the evolution and stability of 
the fluctuations, comparing eventually the decay rates with the Hubble rate to 
demonstrate how a realistic system might actually behave. 
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The high temperature limit of the effective action of the X(p 4 -model has 
recently been considered at one-loop level by Moss et al. [3] (see also [4]). There 
are several reasons for reconsidering their work. First we want an independent 
check of their results. They used an approach based on the local-momentum 
space method for curved spacetime introduced by Bunch and Parker [5], while we 
calculate it in a much simpler way through Feynman diagrams. Secondly, in [3] 
the effective action was expanded in powers of derivatives about a constant field. 
It is known [2] that the finite temperature correction to the two-point function is 
non-analytic at p^ = 0, having different limits when, for instance, {po = 0, p — > 0} 
and {po — > 0, p = 0}. It is, therefore, in general not meaningful to expand about 
the origin, but only in fixed directions in the (po, p)-plane. These expansions 
can be used to analyze e.g. space or time independent solutions. Recalling that 
the variation of the effective action gives an equation of motion for (<p(x)) we 
shall estimate the Fourier spectrum of a solution to the equation of motion. As 
an example we might consider fluctuations about the minimum of the effective 
potential and therefore shift the field by a constant to get to the minimum. Then 
the deviation from the constant field satisfies the on-shell condition p$ — p 2 + m 2 , 
so for a spatially almost constant field we could expand around {po = m, p = 0}. 
Here we have neglected corrections from interactions and quantum fluctuations 
but they have to be small in order for perturbation theory to work. If we were 
interested in other solutions than plane waves (such as solitons) we should study 
the size of the derivatives anew. 

A third reason for studying the A(/> 4 -theory anew is the existence of a 
two-loop diagram of order A 2 which dominates over the one-loop diagram at high 
temperature and when the constant part of the background field, cr c , is close to 
zero. The one-loop tadpole in the X(p 4 -theory, depicted in Fig. 1, goes like AT 2 
for high T but it is momentum independent so that the only one-loop correction 
to the kinetic term comes from the diagram in Fig. 2. It has a high T behavior 
like A 2 cr 2 T. A naive estimate of the T dependence of the two-loop diagram in 
Fig. 3 is A 2 T 2 (though we shall see in Sect. 3.2 that the mass correction goes 
like A 2 T 2 lnT) so it dominates for large enough T. Also, it turns out that an 
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on-shell imaginary part of the effective action, which is important for studying the 
evolution of perturbations, first arises at two loops. 

The paper is organized as follows. In Section 2 we present the derivative 
expansion and the generalized tadpole method for computing the effective action. 
In Section 3 we evaluate the effective action at one and two loop orders for small 
non-constant fluctuations about a constant background. We also comment on 
the existence of a derivative expansion of the effective action. In Section 4 we 
discuss the dispersion relations and the decay rates for the fluctuating modes. We 
show that although most fluctuations will decay during the course of the cosmic 
expansion, some will still be present at the phase transition. Section 5 contains 
our conclusions and a discussion of some open problems. 

2. Methods for computing effective action 
2.1 The derivative expansion 



The effective action r[0] is the Legendre transformation of the generating 
functional W[J] 

e iW[J\ =Af J D$e iS[q?]+i J d4xJ (x)®(x) ^ 

where S[<f>] is the tree level action of the quantum field and M is a normalization 
constant. Thus the effective action of the classical background field is defined by 

m=W[J]- J d*x 6 -^±J(x) , (2) 

where J has been eliminated using the definition of classical field <p(x) = jjj^y, 
assumed to be invertible. From a calculational point of view it is often better to 
express the effective action by an equivalent functional integral form 



= M J £>$ e *s[*+*]+iJV*.7(*)* ; (3) 



where now J(0) = — as can be directly verified from the definition of effective 
action. One can then expand the effective action F[(f)] about an arbitrary field 
a{x): 

oo p n 

r[0(x)] = E - / d ^ r(n) to ^ II W**) - a ^ ' ( 4 ) 



n=0 

where the Taylor coefficients of the expansion are given by 

r<»V;s,]= ^ 



. . .<ty(a;„) 



(5) 



which are fPI Green functions with a background field a. In general r^ n )[a;Xi] is 
a functional of and a non-local function of a^. If we want to write it as 

a local function we need to expand in the derivatives of the fields as well. We 
should however remember that the effective action is genuinely non-local both in 
time and space (e.g. through logarithmic terms from loop corrections). 

Let us assume for the moment that the effective action can be written as a 
derivative expansion of the classical field <f>, which is equivalent to the statement 
that the Green functions are analytical at the origin, when expressed in momentum 
space. Then 



CXJ „ 

r[0(:r)] = J2 d A xY 2k (^ d<j>, □</>, . . .) , (6) 

k=0 J 



where each term T 2k [<j>, 9(f), □</>,.. .] is of the order 2k with respect the 
partial derivative operators. T 2 can always be written in the unique form 
T 2 = \Z\$\d ^fyd* 1 all functionals with two derivatives can always be cast into 
this form by partial integration [4] (note that T 2 cannot in general be written as 
2^(0)0n0). We may thus write down the effective action to second order in 
derivatives by noting that by definition Z and V e ff = —T^°>[(j)] do not depend 
on derivatives, and thus their functional forms do not depend on whether is a 
constant or not. The easiest way to compute them is therefore to consider a 
constant background field a c and afterwards replace a c — > in V e //[a c ] and Z[a c ]. 

It is possible to follow a similar procedure for an arbitrary k (and also 
for fields other than scalars). One must first find a basis, in which it is 
possible to express all the terms containing k derivatives in a unique fashion. 



For example, in the case k = 4, we find that the basis is formed by the set 
{d^d^^d^d^cj), d^d^d^cfrdvcj), d^d^d^d^cj)}. Thus T 4 is uniquely represented 
by the linear combination 

r 4 = r 4 ,i[0]a M 0^0^0a> + r^d^d^d^ + r^d^d^d^ . (7) 

Clearly, different physical problems warrant different expansions. The derivative 
expansion is suitable for e.g. slowly varying background fields (but possibly with 
large amplitudes), whereas modes that approximately obey the on-shell condition 
are best described in terms of an expansion in small perturbations. 



2.2 The tadpole method 



A frequently used method for computing the effective potential is to first 
calculate the tadpole of the theory where the field has been shifted by a 
constant, a c , and then integrate it with respect to a c [6]. In this way an infinite 
class of diagrams is resummed. The reason is that the mass gets a contribution 
from shifting the field and that contribution can be resummed to all orders using 
the effective cx c -dependent mass. We shall here extend the tadpole method to the 
whole effective action. 

Suppose the effective action is expanded about some field a(x) (we consider 
only a real scalar field for simplicity, but the method can easily be extended to 
other fields). Since r[0(x)], given by Eq. (4), is independent of the expansion 
point a(x) we can take the functional derivative of Eq. (4) with respect to cr(y) to 
obtain 



(8) 

T^[a(x);x t }J2^j ~ v) J[(<Kxi) ~ = • 

If we put 4>(y) = a(y) in Eq. (8) and use r[cr(x)] = r^°^[cr(x)] we get 

^™ =!<»[<,(,);,], (9) 



which is the tadpole equation for the effective action. The right hand side of 
Eq. (9) is just the tadpole for a theory where the field has been shifted by a 
non-constant field a(x). 

To see how the calculations work out in practise let us consider the Lagrangian 

£=I(^W-mV)-^0 4 • (10) 

After shifting the field to (p(x) + a(x) we get effectively a space-time dependent 
mass term that cannot be dealt with exactly. This means also that we cannot 
resum the same infinite set of diagrams as for the effective potential. To proceed 
we have to consider a perturbation series in the non-constant part of a(x). Thus 
we write 

a(x) = a c + b(x), (11) 

where a c is constant and b(x) is small. Higher orders are not important 
when the fluctuation energy is dominated by the non-interacting part or when 
b 2 ^ (p 2 +m 2 )/A, where at finite temperature one should replace m by the plasma 
mass m(T). (When there is no risk of confusion, we shall denote the Fourier 
transform of b(x) simply by b(p)). 

The Lagrangian for the shifted field is 

C(cf>(x) + a(x)) =C(a(x)) + ^y<K*) + \ (W*)) 2 - (™ 2 + ^)0 2 (*)) 

,/ n b 2 (x). ,o. . X(a c + b(x)) ,o. . A , 
- ^Mx) + -f 1 )^^) - y c s] y V (*) - , 

(12) 

where r c ^[o"(a;)] is the classical action for the Langrangian in Eq. (10). Because 
of the space dependence in b(x) there are some new Feynman rules that do not 
conserve momentum. The rules are given in Fig. 4. 

Note that the Fourier transform (b 2 )(q) is not equal to (b(q)) 2 but is the 
Fourier transform of b(x) 2 . If we can choose J <i 4 x6(x) = 0, then b(x) only 
includes symmetric fluctuations around a constant field. Then we still have 
(b 2 )(q = 0) = J d 4l xb 2 {x) > if b(x) is not identically zero but measures the average 
size of the fluctuation. In a diagram with all the external momenta equal to zero 
(i.e. for the effective potential) there still remains a contribution from (b 2 )(q = 0). 
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Let us now take a look at the structure of the perturbative expansion of the 
effective action. We can consider an expansion in A, b and h, but since we are 
interested in quantum effects we set the loop expansion parameter h equal to one. 
Higher loops are also suppressed by higher powers of A. The tree-level tadpole 
just gives the classical action. At one- loop and to 0(b) we have the diagrams in 
Fig. 5. The first diagram can be integrated with respect to a c and we get the 
usual one- loop effective potential V e ff(a c ). If we expand the vertex function in p 
in the second and third diagrams we have to zeroth order in p and first order in b 



t,(1) r / \ i „(1) r i f d 4 k (\ A 1 AV2 \ 

<UooM*)\-<Uoo^\ - J (2^)4 [ik^W + 2(Mf ) { ~ P) 



(13) 



Id f d A k \a c . . . . d 2 Tr , ... , 



2 8a c J (2n) 4 k 2 -M 2 
where 

M 2 = m 2 + Xa 2 /2 . (14) 

This is what we get if we first replace a c by a c + b(x) in dV e ff(a c )/da c and then 
expand to first order in b. We can easily understand this from the derivative 
expansion of the effective action, as was discussed in Sect. 2.1. 

The tadpole diagrams that are linear in b(p) are, of course, essentially the 
two-point functions for the theory with shifted a c . In the following we calculate 
instead 1?( 2 ) directly and one can construct the effective action from Eq. (4). The 
effective action in momentum space to 0(b 2 ) is 

T[a^b] = -V eff (a c )+TW(a c )b(k = 0) + l J -j0L b (-k)r™(a c ; k)b(k) , (15) 

and rW is zero at the minimum of V e //. 



3. Finite temperature effective action 



There are at least two formalisms for calculating vertex functions at finite 
temperature, the imaginary and real time formalisms (ITF and RTF). They should 
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both give the same result when used correctly. In ITF there is a naive way of 
extracting the leading high T behavior by only including the mode with zero 
Matsubara frequency (n=0) (for bosons). For the diagrams under consideration 
here this method is not reliable. The correct method should include a summation 
over all Matsubara frequencies and an analytic continuation at the end. RTF is 
easy to use for one-loop diagrams but for the two-loop diagrams we have employed 
instead ITF with the convenient method developed in [7-9]. 



3.1 One-loop diagrams 



Let us start by giving the well known result for the zero temperature part. The 
effective potential can be found in [10] and we repeat it here for completeness: 



VeffM = fm 2 + + ^ In ^— ) + .4 + Ba, + ( o; 
The renormalized two-point function is (see Fig. 2 



r?-W*(*)] = j^f 



2^2 



X 2 a 



AM 2 — p 2 
-2\l ^ arctan 



AM 2 - p 2 



In 



M 2 
Att^j? 



(16) 



(17) 



AV C 2 

32tT 2 

AV C 2 

32tT 2 



-2 + 



P 2 


-In ^ 


6M 2 


— )- 


2* 
l 3V3 



M 2 
47^ 



+ 0(p 4 ) , 



J9 2 - M 2 

} M 2 



-In 



M 2 

471"^ 



+ £>((p 2 -M 2 ) 2 ) 



While the expansion around p 2 = agrees with [5,11], we should like to point out 
that the expansion around p 2 — M 2 is actually better suited for solutions close to 
the mass shell. 

The finite temperature effective potential has been calculated in several papers 
(see e.g. [12]) and we just give the result (see Appendix A for the definition of Ff) 



V eff A°c) = -^Ft{T,M 2 ) 



T 4 tt 2 T 2 m 2 a 2 XT 2 TM 3 (a c ) 



90 



+ 



24 



+ 



2 24 



12tt 



(18) 



The dominant mass correction at high T is XT 2 / 24. 

Next we look at the effective action at finite T and here we notice some 
problems if we estimate the high T limit using only the zero Matsubara frequency 
mode. In that approximation we would find 



r ?,Lo(Po,p) = 



X 2 a r T 1 



16/T p 



arcsm 



P 2 +Po 



— arcsin 



pV4M 2 + p 2 -2p2 +p 4 /p2 



P -Po 



(19) 



pV4M 2 + p 2 -2p2 +p 4 /p 2 



where p = |p|. This function has a well defined limit when p^ — > and a threshold 
at {po = M, p = 0}. These properties are however not shared by the correct 
two-point function. In particular the threshold is unphysical since there is no 
decay process allowed at that momentum [13]. To get the correct answer we can 
use either ITF or RTF to obtain 

f)f^ ^_ AV c f°° dk 1 lT1 ( (go - P 2 + 2pfc) 2 - 4p 2 c^ 

where 



^ (P0 ' P) -"l6^y w F ln V(^-P 2 -2pfc) 2 -4p 2 u; 2 ), ) ' (20) 



If we set po = 0, expand in small |p| and take the high T limit we get 



r^ ( o,p) = -^|r^M?in 



o w p 



p - 2/c 



p + 2/c 



\2 n .2 T-i ^2 
c M _ P ^ 



(21) 



16tt M v 12M 2 ' 

where the p-term coincides with the result of [3].* A possibility of avoiding the 
problem of non-analyticity at the origin is to expand around the mass shell. When 
we expand around {pi = M 2 , p = 0} we find 




Note that there is a misprint in Eq. (24) in [3]. 



One can, of course, obtain higher derivative terms by simply expanding to higher 
order in (pjj - M 2 ) and p 2 . 

The expansion in Eq. (22) could be used to write down the dispersion relation 
near mass shell. We shall return to dispersion relations in Sect. 4. After 
performing the resummation of the tadpole to get the thermal mass of AT 2 /24, no 
particularly interesting feature is found in the high T limit. The only appreciable 
effect is the thermal mass. For instance, on-shell there is no imaginary part at 
the one-loop level. Such an imaginary part, which is of essential importance when 
one wants to study the evolution of the perturbations, is first generated at the 
two-loop level. 

3.2 Two-loop diagrams 

For T/M large, higher loop diagrams go as higher powers of T than the 
one-loop diagram. On the other hand they are also suppressed by higher powers 
of A. At two-loop order the diagrams in Figs. 3 and 6 are of the order A 2 , all 
other are of higher order. The leading diagram in T/M to each order in A is the 
one with as many tadpole insertions as possible. As can be seen from Eq. (18) 
it gives a AT 2 /24 contribution which can be resummed by just shifting to an 
effective temperature dependent mass M 2 + XT 2 / 24. This is well known at the 
one-loop level [12] but at the two-loop level one must be careful not to double 
count diagrams. Parts of the leading "double-bubble" (Fig. 6) are already included 
in the one loop tadpole when the mass is shifted. An easy way of keeping track of 
the counting is to introduce a finite, T dependent counter term as in [14]. In the 
"setting sun" (Fig. 3), however, the leading corrections can be resummed using a 
T-dependent mass. In this Section we perform the calculation without resumming 
tadpoles, and defer the resummation to Sect. 4. 

The diagram in Fig. 6 does not depend on the momentum so it can be 
absorbed in the effective potential. It has been calculated both at zero [11] and 
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finite [15] temperature. As for the one- loop case we state the result for the 
effective potential for completeness 




+ A + Ba 2 c + Ca\ ■ 



At finite T the result is 



(24) 

Here v'l^o^) is a second order polynomial of a c that depends on the renormalization 
condition at one-loop, and F?(T, M) and I(M/T) are defined in Appendix A. We 
have not verified the expressions in Eqs. (23) and (24), which can be found in 
[11,15], but instead concentrated on the high temperature expansion of Eq. (24). 
Most of it was derived in [15] but here we also show that (see Appendix A) 

T 2 ^TT 2 /M 2 \ 

/(r - M) ^?r ln (^) + 0(r2) ' (26) 

There is also a subleading term that goes like T 2 which need not be numerically 
small compared to the one given in Eq. (25). 

When it comes to the effective action the finite temperature part of the "setting 
sun" is particularly interesting since it dominates over the one-loop "bubble" at 
high temperature T 3> cr 2 /M. In the "setting sun" there is a T dependent UV 
divergent constant term but that cancels against the "double-bubble" diagram so 
that there are no T dependent infinities when all 0(X 2 ) diagrams are included. 

The derivative term in the "setting sun" is finite (since there is no p-dependent 
term in the "double-bubble" that can cancel infinities) and one might be tempted 
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to try high T expansion using the zero mode approximation. Expanding in po and 
p one finds, after some algebra, that 

(2) _ A 2 T 2 f d 3 k dH 1 1 1 



6 J (27r) 3 (27r) 3 k 2 + M 2 (k + l) 2 +M 2 -p 2 + (p-l) 2 + M 2 
A 2 T 2 / . , 1 ^ (26) 



576tt 2 M 2 \ A + p o- 9^ 

where A is an UV divergent momentum and temperature independent constant. 
Our experience from the one-loop calculation in Sect 3.1 shows, however, that we 
cannot trust the zero mode approximation. In fact, although the imaginary part 
of 1?( 2 ) at {po = M + ie, p = 0} is UV finite, when it is calculated from Eq. (26) it 
is zero, whereas the correct value is [14] 

Imr( 2 )(p = M, P = 0) = ^ . (27) 

Also, the leading T dependence of the "setting sun" goes like T 2 lnT as we shall 
show below. We must therefore sum over all Matsubara frequencies and the easiest 
method to do that is the one described in [7-9] and [16]. 

The two-point function can then be written as 

r^ oo >0, P ) = y / ( 2^3 (2^3 jf dre^A(-r, k) A(-r, 1) A(r, p + k + 1) , (28) 
where 

1 e^ k i 1 

A(r, k) = — - (-f - e -"*M + -e" k l r l) . (29) 

v ' 2u; k - 1 - 1 ' y ' 

After carrying out the r-integral we find 



r (2) 2 -Zoo P (po,p) = -y j 



d 3 k dH 1 1 



) 3 (2tt) 3 u k ui pi - (u; k + uj\ + w p+k+ i) 2 

A^ r d 3 k d 3 i l ^ l 

T J (2^(2^^ /s(u;k) \ (p ±^ k ) 2 -( Wl + Wp+k+1 ) 2 (3°) 

A 2 f d 3 k d 3 l 1 w , ^ 1 

~o~ / To ^3/o ^3 Jb\^\l)Jb\^\) > 7 — 77 77 — 75 2 > 

8 J (2tt) 3 (2tt) 3 uj^ujx ^ (p ± cu k ± u)\Y - w^ +k+1 
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where the summations are over the two and four possible combinations of ±. We 
keep only the terms with one and two distribution functions /s(cu) since we are 
only interested in the finite temperature part. The term with one distribution 
function has a temperature dependent, but momentum independent, infinite part 
which cancels when all contributions are added. We also subtract this part to get 
a finite answer. All the angular integrations can then be performed explicitly but 
we leave one integral for notational simplicity. ^2-ioop(Po^ p) nas a l so an imaginary 
part which we determine through an analytical continuation corresponding to the 
retarded two-point vertex (p — > p + ie). The final result is 



r 2-Zoo P (P0,p) =- 



A 2 



dk dl kl 



f B M(G?\po,p) + 5G? ) (p ,p)) 



8(2tt)- J w k wi 
+ f B (^)fB(ooi)G^(p ,p) 



(2), 



(31) 



where 5G\ (po,p) regulates the infinity [14], and we have separated the terms 
with one and two distribution functions: 



pk+p 




/ *E< 




l\k-p\ ± 





pI - (±w k + uj\ + u>i- z y 



G { 2 } (po,p) 



pi - (±w k + u\ + ui +z ) 2 
+ in [9 (pi - (±u; k + ui + uji- z ) 2 ) -6(pl- (±w k + ui + uji +z ) 2 )] }> , 

(po ± w k ± co\) 2 - m 2 -(z- I) 2 



±- l k+v dzy[\n 

2 P J\k-p\ t%\ 



(po ± cu k ± cui) 2 - m 2 - (z + l) : 



+ in sign(p ± <^ k ± u\) [9 [{z - I) 2 + m 2 - (p ± cu k ± cui) 2 ) 
-9((z + l) 2 + m 2 -(p ±iv ]i ±iv l ) 2 )]^ , 

_( 2 ), , 2/ [ k +*dz 
$G\ ; (p ,p) = - / — • 

P J\k-p\ z 

(32) 

It is interesting to note that the two-loop result Eq. (31) is analytic at the 
origin, unlike the one-loop action Eq. (20). The leading mass correction has 
been calculated explicitly by Parwani [14] at high T and at {po = M, p = 0} (see 
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Appendix A for details): 

r& oo > = m, p = o) ^ -\-^ 2 t 2 Hy) ■ ( 33 ) 

One can now construct the effective action using Eq. (15). 

4. Evolution of perturbations 

When studying the evolution of perturbations about a constant field, the 
knowledge about the two-point function is usually not enough. The action has, 
apart from some special cases, also a linear part, so that the equation of motion is 
of the form T^ 2 \a c )b = — T^(a c ). The right hand side of this equation can be 
eliminated by writing the field in the form b(x) = bh(t) +Sb(x). The first part is 
only time dependent and satisfies the equation of motion with vanishing spatial 
derivatives. It describes the homogeneous rolling of the field on the effective 
potential surface as long as bh(t) remains small. The second part, 5b, describes 
the fluctuations on the homogeneous field, and satisfies the dispersion relation. 

As an application we shall now consider the evolution of local perturbations 
b(x) around the high T minimum a c = 0, which is the initial state for the phase 
transition, and for which = 0. The one-loop contribution to , Eq. (20), 
vanishes, and the lowest order correction is provided by the two-loop diagram of 
Fig. 3 as given by Eq. (31). For small b(x), when we may neglect terms higher 
than 0(b 2 ) in the expansion for the effective action, the evolution of different 
Fourier modes is independent and given by the dispersion relation 

Pi = P 2 + m% + X R (p , p; to) + iE 7 (po, p; to) , (34) 

with 

E fl (p , p; m) = - Re (t™ (po, p; to) - (to, 0; to)) , 
£j(po,p;to) = -Imr (2) (p ,p;TO) , 
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(35) 



where the renormalized mass tur has been introduced to absorb the order C(A 2 ) 
mass correction, part of which comes from (m,0;m). We have also indicated 
explicitly the mass-dependence of the two-point function. We may write 

AT 2 

m 2 R = m 2 + — + 0{\ 2 ) ee m(T) 2 + 0{\ 2 ) , (36) 

where m(T) is the familiar one-loop plasma mass. In fact, to order 0(X 2 ) we may 
simply replace the mass parameter m by m(T) in 1?( 2 ) (p , P; m). This amounts 
to plasma resummation of the propagators in the loop diagrams. (As discussed 
before, there is no double counting in this case). Thus to lowest order we may 
write, denoting Impo = — 7/2 and Repo = u, the dispersion relation as 

^ 2 =p 2 + m(T) 2 -i 7 2 + E fl ( W ,p;m(T)) , 
7 = - £/(u;,p; m{T))/u . 

We may solve Eq. (37) approximately by setting uj 2 ~ cu 2 = p 2 +m(T) 2 . To this 
end we have studied the behavior of E^(p ,p;w) and E/(p ,p;w) numerically. It 
turns out that when p — > 00, the contribution from the part with two distribution 
functions becomes unimportant. Roughly, numerically we find that 7 ~ T 3 /p 2 
which is also an upper limit that can be derived analytically. We also find that for 
large p the real part is independent of p and goes roughly like (27r) _4 A 2 InT, so 
that for A 1 this part is never important for the dispersion relations. The limit 
p — > 00 is thus completely dominated by the part with one distribution function, 
+6G1 in Eq. (31). Scaling the relevant part of Eq. (31) by p and taking the 
limit, it is possible to show that (see Appendix B for details) 



r (2) K,p)^ g(£p Fl2(T ' m(T)) 



1 p2 4- ■ ' 



(38) 



where F 2 has been defined in Appendix A, and for large T and A< 1, we have 
F 2 ~ 7r 2 T 2 /6. In the high temperature limit the dispersion relation Eq. (37) 
for the real part reads then pi = p 2 + m(T) 2 + 3(2yr)- 4 Am(T) 2 ln(p 2 /m(T) 2 ) so 
that T,it(po,P',m) never contributes significantly. Thus u = uj p is a self-consistent 
approximate solution. 
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As discussed in [13], 7 is the rate by which a given plane wave mode 
thermalizes. We have computed Im T^ 2 ' from Eq. (31) numerically, and the 
resulting thermalization rate is drawn in Fig. 7 for a few temperatures. The 
imaginary part exists for all four-momenta obeying the (approximate) dispersion 
relation, and this means that every small perturbation b(x) about the high T 
minimum a c = will eventually smoothen out and thermalize. 

The essential question then is, how rapid is the rate of thermalization as 
compared with the expansion rate of the Universe, given by the Hubble parameter 
H. Assuming the Standard Model degrees of freedom, at high temperature we 
may write 

\ 90 J Up- M P • ( ' 

where M P = 1.22 x 10 22 GeV is the Planck mass. Comparing (38) and (39) we find 
that all momenta obeying 

p < 2.8A 2 x 10 14 GeV , (40) 

decay in less than a Hubble time. 

Recall that the critical temperature and the intrinsic mass scale (the inverse of 
the wall thickness of a bubble) are typically related by V\T C ~ M. This means 
that perturbation modes of the size of the critical bubble have not yet decayed 
by the onset of phase transition provided A < 10~ 5 (M/100 GeV). For such 
small A the inherent local perturbations would certainly be important for bubble 
nucleation. 



5. Conclusions 



We have found that the finite temperature effective action can be relatively 
easily computed using the tadpole method generalized for space-time dependent 
fields in Sect. 2.2. In that case resummation is however not possible except for 
constant background fields a c . Nevertheless, one can do perturbation theory in 
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small fluctuations b(x) about a constant background, which in effect amounts to 
computing Green functions for the shifted theory. To lowest non-trivial order one 
then just evaluates the two-point function, which we did to two-loop order in 
Sect. 3. After that it is straightforward to write down dispersion relations and 
to study the evolution of the fluctuations. Since this only requires an expansion 
of the effective action about the mass shell we do not encounter the problem of 
non-analyticity at p M = 0. We also observe that the zero Matsubara frequency 
approximation for the leading high temperature terms is not correct in general. In 
our case, even the UV-finite imaginary part needs a summation over all Matsubara 
frequencies. 

Our main conclusion is that small space-time dependent perturbations about 
the high temperature minimum do not automatically thermalize before the phase 
transition starts. This would naturally affect bubble formation in a first order 
phase transition, although to what extent is not clear. In spite of the fact that 
this conclusion holds strictly speaking only for the A^ 4 -model, we could view the 
A</> 4 -model as an effective theory of the order parameter of a more realistic theory, 
such as the Standard Model. The gross features of the present study would then 
presumably show up also in the full Standard Model. 

It is however clear that the estimate at the end of Sect. 4 is very rough and has 
to be improved for an application to a real physical system. The main issue is the 
strength of the effective coupling. Some recent estimates seem to indicate that the 
temperature dependence can be quite strong, and that the four-point interaction 
for zero momentum particles goes in fact to zero at the critical temperature of 
a second order phase transition [17,18]. This is caused by the vanishing mass at 
the critical temperature. The coupling constant decreases at the same time as 
the mass and thus the relevant p also decreases. It is not clear what the net 
effect of the coupling constant renormalization is. Furthermore, the electroweak 
phase transition is likely to be first order, albeit possibly very weakly so. These 
questions can only be addressed after a more detailed study of the damping rates 
in the full electroweak theory [19]. 

We have only considered small perturbations. It is conceivable that large 
perturbations could, because of the non-linear coupling of the modes, contain 
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instabilities so that some perturbation modes would actually get amplified. Further 
studies in this direction are certainly needed. 



18 



Appendix A 



In this Appendix we give some useful expansions employed in Sect. 3 and 
calculate the leading T dependence of Eq. (31) as well as the function I(T/M) in 
Eq. (24). 

The function 

r°° Ah Um 

iC(T,M) = / — ^(a/) , u = aA 2 + M 2 , (Al) 
Jo u 

often occurs in loop calculation results at finite temperature (cf Eqs. (24), (38)). 
It satisfies the relation 

^ = -^K-\ ™>1. (A2) 
In this paper we need F 2 and Ff and they have the high T expansions 



1V ' ; 6 2 4 y AnT ! 

Ft(T, m) = tEL - friMl + ^IM1 + IM1 M EL ) + 

1V ' ; 15 4 2 16 y AnT ! 



(A3) 



Let us now turn to the leading T dependence of Eq. (31) and the function 
I(T/M) in Eq. (24). The leading high T expansion of Eq. (31) was derived in 
[14] . We have checked the coefficients of the T 2 In T terms in the part with two 
distribution functions, using a different method, and found them to be correct. 

The essential integral in the part with two distribution functions of Eq. (31) is 



T*h(T/M) = ^ y^f B fa)f B {un) In 
Jo ^k^l 



k-l 



k + l 



{AA) 



where a = M/T and {po = M, p = 0}. We expect I\ to go like In a for small a so 
we compute the derivative with respect to a 2 . Using 

d f f(vx)f(v y ) \_ f(u>y) 1 d f f(u x ) \ | fjux) 1 d f f(u y ) \ 
da 2 \ uj x ojy J uj y 2x dx \ uj x ) uj x 2y dy \ uj y ) 
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where uj x = Vx 2 + a 2 and f{u) x ) = (e u>x — 1) 1 , we find after partial integration 

dh f°° dxdy ... . In 2 



f dxdy 1 

/ /K)/K ^ — 

Jo u x uj y a 2 



da 2 ] n LUx.ojq, a 2 4 

for small a. After integration with respect to a 2 we get 

2 rp 

T 2 h(T/M) = — — T 2 ln( — ) + C(T 2 ) . (A6) 

£ 1VJ. 

The function I(T/M) in Eq. (24) can be expanded in a similar manner. It is 
defined as [15] 



I(T/M) = %L r ^y fM f(u y ) In 

M Jo U x UJ y 

The same trick works and we find 

5 7i~ 2 T 2 . T 



(a 2 + 2xy) 2 - Auj 2 x uj 2 y 



(a 2 — 2xy) 2 — 4:U 2 u) 2 



(AT) 



/(r/M) = -- T ^ln(^) + 0(^), (AS) 

where the factor of 5/6 comes from 

r°° dx dy 2x 4 - x 2 y 2 +2y 4 _ 5 n 2 

J (x 2 + l)(y 2 + l) (x 2 -xy + y 2 )(x 2 + xy + y 2 ) ~ 6 2 

It is the same diagram, the "setting sun" , that leads to I and ii but I is computed 

at po = since it occurs in the effective potential while ii is computed on-shell, 

Po = M. We therefore do not expect them to be equal but only of the same order 

of magnitude and to have the same sign. 



Appendix B 

The asymptotic behavior of 1?( 2 ) can be calculated from Eq. (31) for all 
temperatures in terms of the functions F™. In the effective action the part with 
one distribution function, + 5G^ in Eq. (31), dominates over and thus 

the behavior of the two-point function in the limit p — > oo is solely determined by 
Gf + 5G ( ? ) . We define 

Jo J\k-p\ 2 P ± l^i Po + K ~ ( ±w k + wi + c^ +2 ) 2 WiJ 

(SI) 
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Ii can be evaluated in the limit p — > oo. After scaling the integration variables 
k, I and z by p we can extract the leading behavior of the integral 1^. The 
logarithmic piece has now a finite, non-zero limiting value, and to the leading 
order we can take this limit in the integrand. We obtain 



f°° dkdlk 3 
Jo U p k 

/■ 1+fe ,lv-f n l + ie-(l + \l-z\±k)' 2 4z 

x / dz-y < In —— - ' H ; 

2^1 l + ze-(2/ + ^±/c) 2 + M 2/ pS 



(52) 



Note that we cannot take the limit in the whole integral, because both the k and 
/-integrations are singular at the origin when p — > oo. The easiest way to proceed 
is to first perform the / -integration and after that the ^-integration. This yields 



/•oo 3i, 



p 2 

2/cln-^r + i2nk 
M 2 



which, apart from prefactors, is 1?( 2 ) in Eq. (31). By rescaling k — > k/p we obtain 
the result Eq. (38). The logarithmic behavior of the real part of the two point 
function is due to the IR-singularity in the UV-regulating term. 
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Figure Captions 

Fig. 1 The one-loop tadpole. 

Fig. 2 One-loop diagram for T^ 2 \ 

Fig. 3 Two-loop diagram for 1?( 2 ) ("setting sun"). 

Fig. 4 Feynman rules for momentum-dependent external legs. 

Fig. 5 One-loop diagrams for the effective action in the generalized tadpole 
method. 

Fig. 6 Two-loop tadpole ("double bubble"). 

Fig. 7 The thermalization rate 7 for on-shell configurations uj ~ uj v for 
T= 10 m(T) (dash-dotted curve), 100 m(T) (dotted curve) and 1000 m(T) (solid 
curve) . 
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Abstract 

The effective action is computed for the A</> 4 -theory at finite temperature for 
small perturbations about a constant background field, using a generalized tadpole 
method. We find the complete effective action, including the real and imaginary 
parts, to all orders in derivatives and to order 0(X 2 ). We demonstrate that the 
high T approximation, where only the zero Matsubara frequency is included, is 
incorrect for the imaginary part even though it is UV-finite. The solutions of the 
dispersion relations show that initial perturbations do not necessarily thermalize 
fast enough to be absent at the onset of phase transition. 
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